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1 Introduction 

Let (r, cr) be the spherical coordinates in M.^ . If n is a harmonic function in \ {0} 
written under the separable form 

u{x) =r-^Lj{a) (1.1) 

it is straightforward to check that uj is an eigenfunction of the Laplace-Beltrami operator 
—A^j^_-^ on the unit sphere S'^"^ C and /? is a root of 

- (iV - 2)X - A = 0, (1.2) 

where A > is the corresponding eigenvalue. The function uj is called a spherical harmonic 
and its properties are well-known, since such functions are the restrictions to the sphere of 
homogeneous harmonic polynomials. More generally, if Cs C is the cone with vertex 
and opening S C S'^~^, there exist positive harmonic functions u in Cs under the form 
(jl.ip which vanish on dCs \ {0} if and only if /? is a root of (jl.2p . where, in that case, 
A := Ag is the first eigenvalue of — A^^^ -^ in Wq''^{S). These separable harmonic functions 
play a fundamental role in the discription of isolated interior or boundary singularities of 
solutions of second order linear elliptic equations. If the Laplace equation is replaced by 
the j?-Laplace equation 



Apu := -div {\Du\p-'^Du) = 0, (1.3) 
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{p > 1), the same question of existence of separable p-harmonic functions, i.e. solutions 
of (|1.3p in the form (jl.ip . was considered by Krol [llj, Tolksdorf ^19j, Kichenassamy and 
Veron [lOj. If u in (II. Ih is p-harmonic, then the function co must be a solution of the 
spherical p-harmonic equation, 

- div (^{P^u^ + |V'a;p)P/2-iv'a;) = - 1) +p - TV)!/?^^^ + \Vuj\y/^-^uj, (1.4) 

on S"^^^, where V and div are respectively the covariant derivative identified with the 
"tangential gradient" and the divergence operator acting on vector fields on S"^"^. Two 
special cases arise when either p = 2 or = 2: if p = 2, (jl.4p is just an eigenvalue problem 

-A'uj = P{P + 2-N)lu, (1.5) 
where A' is the Laplace-Beltrami operator on S^'^. When N = 2, equation ()1.4p becomes 
- (^{PW + \LOe\y/^-'ioe)^ = PiPip - 1) + p - 2KPW + |u;e|2)P/2-i^, (1.6) 

where 6 S [0, vr]. Introducing the new unknown cp := LOg/uj, (II. 6p is transformed into a 
separable equation, 

- ((/?' + <pY'-'<p)^ = Up - 1)0' + PiPip - 1) + P - 2)) iP^ + 02)^'/2-i. (1.7) 

This equation was completely integrated by Krol [H] in the case P < 0, and Kichenassamy 
and Veron [10] in the case P > 0. It turns out that for any integer k > there exist two 
couples {Pkjfpk) and {Pk,4>k) where Pk < 0, Pk > 0, and (j)k and (pk are anti-periodic solu- 
tions of the corresponding equation (|1.7p . Furthermore (pk and (/>fc are uniquely determined, 
up to an homothety. 

A remarkable breakthrough was realized by Tolksdorf |19j when he proved that for 
any smooth domain S C S'^'^ there exists a couple (/?,(/>) where P < and (p G C^{S) 
is positive in S, vanishes on dS and solves (II. 4p in 5". Furthermore P := Pg is unique 
and is determined up to a multiplicative constant. Tolksdorf's result is obtained by 
constructing a p-harmonic function u in the cone Cg generated by S with a compactly 
supported boundary data and by proving, thanks to a kind of Harnack inequality up to the 
boundary, the "equivalence principle" , that the asymptotic behaviour of u is self-similar. 
Later on the existence of a couple {P,(p), with P := P^ > and (p, as above, positive 
solution of (11. 4p in S vanishing on dS is proved by the same method in [2Tj, therefore we 
shall refer to the two cases P > and P < as Tolksdorf's results. The structure of these 
spherical p-harmonic functions is studied in [5] when p = N. These regular (/? < 0) and 
singular (/? > 0) separable p-harmonic functions play a fundamental role in describing the 
behaviour of solutions of quasilinear equations near a regular or singular boundary point 

In this article, we give a new proof of Tolksdorf's results, entirely different from his. 
Actually, performing a change of variable, we embed our problem into a much wider class 



3 



of quasilinear equations. Indeed, if u; G Wq'^{S) is a positive solution of (II. 4p in 5 C S'^ ^, 
which vanishes on dS, then the function v defined by 

V = ——InuJ 

solves 

' -div ((1 + \V'v\^Y^^~^ V'v^ + (3{p - 1) (1 + \V'v\^Y^^-^ |V't;|2 

= - {f3{p -l)+p-N){l + \V'v\Y^^'^ in S (1.8) 
lim^^dS v{a) = oo. 

Notice that this equation is never degenerate and v is (actually C°°) in S and satisfies 
the equation and the boundary condition in classical sense. Our construction of solutions 
of ()1.4p relies on a careful study of the quasilinear problem (jl.Sp , and on the interpreta- 
tion of the constant in the right hand side of (II. 8p as an "ergodic constant" . Furthermore, 
having an intrinsic independent interest, this study will be performed on any compact 
smooth subdomain of a Riemannian manifold, without refering to the p-Laplace equation 
(jl.Sp . Our main result is the following: 



Theorem A. Let {M,g) be a d- dimensional Riemannian manifold with nonnegative Ricci 
curvature, and let V and divg he respectively the covariant derivative and the divergence 
operator on M . Then for any compact smooth subdomain S C M and any P > there 
exists a unique positive constant \p such that the problem 

' -div, ((1 + \Vv\^Y'''~^ Vt;) + /3(p - 1) (1 + |Vt;|2)^/'-^ \Vv\^ 

= _A^(l + |V^;l2)P/2-i ^n S (1-9) 

\va\^^QSv(x) = CQ. 

admits a solution v € C'^(S). Furthermore, v is unique up to an additive constant. 



The result of Theorem A is the typical statement of an ergodic problem, indeed the 
constant can be seen as the unique ergodic constant for the equation obtained after 

dividing by (l + |Vvp)^''^ ^ (see (12. ip ). Observe also that (II. 9j) may be reformulated if 
we set uj = e~^^, then w is a solution of 

r -diVg {{p^u^ + |Va;|2)P/2-iVa;) = pXpiP^uj^ + \Vuj\'^)p^^-^uj in S 
[ UJ = on dS 

When p = 2, problem (jl.lOp reduces to an eigenvalue problem since pXf^ = Aj(S'), the 
principal eigenvalue of the Laplace-Beltrami operator in S. In that case the connection 
between (jl.9p and (jl.lOp dates back to the stochastic interpretation of principal eigen- 
values (see e.g. [16J, [T7]). In the nonlinear framework with p 7^ 2, by proving that the 
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mapping /? i— > A/j is continuous, decreasing and tends to oo as /? — > 0"'", we conclude that 
the equation = {j3{p — l)+p — d—l) has a unique positive solution. As a consequence 
we generalize Tolksdorf's result as follows. 

Theorem B. Under the assumptions of Theorem A, for any compact smooth subdomain 
S of M there exists a unique /S := (3^ > such that the problem 

j -divg ((/32u;2 ^ I Vu;|2)p/2-i Vcj) = /? {p{p - I) + p - d - 1) {p'^oj'^ + | Vcu|2)p/2-1u; in S 

u! = on dS, 

(1.11) 

admits a positive solution uj € C'^{S) n C^(5). Furthermore ui is unique up to an homoth- 
ethy. 

Of course, we obtain similarly that for /3 < there exists a unique (3 := (3g < Q such 
that Xp = {f3{p — 1) + p — d — 1). Tolksdorf's results then follow as a particular case by 
taking {M,g) = {S'^^^ , go), where S^^^ is equipped with the standard metric go induced 
by the Euclidean structure in M.^ . 



2 The singular case 

In the following, we consider a general geometric setting and we recall some elements 
of Riemannian geometry (see e.g. [13], [E]). Let {M,g) be a complete d-dimensional 
Riemannian manifold with metric tensor g = (gij), inverse g~^ = (g^^) and determinant 
\g\. If X and Y are two tangent vector fields to M, we denote by 

X.Y = Y^gijix)X'Y^ 

ij 

their scalar product in the tangent space T^M. Let Xj, j = l,...,d, be a local system 
of coordinates: if n G C^(M), the gradient of u, quoted by Vu, is the vector field with 
components (Vu)* = J2k Therefore 

Vu.Vu = |Vup = ^g'^{x)u^^u^^. 
li X = (X*) is a vector field on M, the divergence of X is defined by 

divgX = (y^\x'') . 

Recalling that, in local coordinates, the Christoffel symbols are 

2 Y V dxi dxj dxj^ ' 
the second covariant derivatives of a function u are 

VijU = UxiXj — ^ij^Xk, 
k 
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while the Hessian is the 2-tensor D^u = (Vyu). Finally, AgU = trace{D^u) = divgVu is 
the Laplace-Beltrami operator on M, locally expressed by 

We denote by Riccg the Ricci curvature tensor of the metric g. In particular, if (M, g) = 
{S^~^,go), then fficc^,, = {N - l)go. 

In all the sequel p > 1 is a real number. We prove next the result of Theorem A, which 
we restate here for the reader's convenience. 

Theorem 2.1 Let S d M he a smooth hounded open domain of M such that RicCg > 
on S. Then for any /3 > there exists a unique > such that there exists a function 
V G C'^{S) satisfying 

-Agv - (p - 2) ^ ^ ^^^^^ + I3{p - 1)|V^|^ = -\p inS ^2.1) 

\\TCLx-^QSv{x) = OO. 

Furthermore, v is unique up to an additive constant. 

Proof. As in the usual approach to crgodic problems, we start by considering the problem 

-Agv, - (p - 2):^^^^^ + P{p - l)\Vv,\^ + ev, = mS ^2.2) 

\im^^gsv^{x) = 00, 
where e > 0, and then we study the limit when e — > 0. 

Step 1: Construction of super and sub solutions. Since dS is C^, the distance function 
p{x) = dist {x,dS), where the distance is the geodesic distance, is a positive function 
is some relative neighborhood ^Sf^ = {x G M : \p{x)\ < 5} of dS; here p{x) is the signed 
distance, equal to ip{x) according x(zSoyx(^M\S. Then \Vp{x)\ = 1 in 'Ns- We 
extend p outside into a C^(M) function p. Next we consider the function 

u{x) = Hp{x)) - Mop{x) + ^ Vx G S, (2.3) 

where the Mj > are to be chosen later on. Then 

Vt2(x) = --^ (1 + f3Mop{x)) Vp(x), 

Notice that this last identity implies 
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Next 



-A„u -ip- 2) ■ ^ = -A„n - — -— — , ,„ 

- I^^l' I ^3~P(^ I .^.V/.^l ^^"^^ \V~p\Hi + m,~P? 

f3~P^ rp rp^ l3^p^ + \Vp\^{l + l3Nhpf 

^ (p - 2) (1 + /? Mo /5) Vp • V [|Vp|2(l + /3Mo p)^] 



2/3/5 /32p2 + |Vp|2(i+/5Mo^)2 

After some lengthy but standard computations, one obtains the following relation 
D'^uVu.Vu 

+ IJ[p - ijlVul" + eu 

(2.4) 



-Agu -{p- 2) + /3(p - l)|V^x|2 + eu 

1 + Vnr 



;51n(p) + 2{p - l)Mo|Vp|2 + ^l;p{x) + Mi 



where ipp is a function depending on /3 (and on Mq), but which remains bounded on 
S, uniformly when /? remains in a compact subset of (0,cxd). Since |V/5| = 1 near the 
boundary, it is possible to choose Mq and Mi such that u defined by (12. 3p is a supersolution 
for (j2.2p . Moreover, Mq and Mi can be chosen independent of /3 whenever it varies on a 
compact subset of (0,oo). 

One finds similarly that the function 

u{x) = ln(p(x)) + Mop(x) - — Vx G 5, (2.5) 



is a subsolution of (j2.2p . with Mg and Mi chosen as for u. Moreover, for < h < 6, we 
can approximate u and u respectively from above and from below by 

Uhix) = ln(p(x) -h)- MoiPix) -h) + (2.6) 

uj^ix) = -i ln(p(x) + h)+ MoiPix) + h)-^, (2.7) 

which are, respectively, a supersolution in {x £ S : pix) > h} and a subsolution in S". 
Together with the comparison principle, these super and sub solutions will be used to 
derive estimates on the solutions of (12.21). 



Step 2: Basic estimates. In this part, by using the classical Bernstein's method ([2j), we 
derive the fundamental gradient estimate for the solutions u G C'^{S) of 



- Agu - (p - 2) "^ ""^^"^9^ + Pip - 1)|V^P + eu = in 5. (2.^ 
1 + Vnr 



We recall the Weitzenbock formula (see e.g. [T]): 
1 



Aj;|Vup = iD^up + V(AgM).Vn + RiccgiVu, Vu), (2.9) 
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and the Cauchy-Schwarz inequality for D^u 



1, 



Let m = inf{i?icCg(Vn, Vn) : |Vn[ = 1} > 0, then 

^ ■ '— l2 \ 1 I A „.|2 , _tv7„.|2 



^AglVnj^ > ^lAgnp + m|Vur + V(Agn).Vn. (2.10) 
If we set z = |Vnp, we can re- write ()2.8p as 

A.^^ = -^^Y^^ + /3(p-l)^ + 6n in 5. (2.11) 
Using the fact that 

V(Vz.Vii).Vii = D^zWu.Vu + ^|Vzp, 

we obtain 

V(A^u) Vu - ~ ^'^Vn.Vtx (p - 2) |Vz|^ ^ (p - 2) (Vz.V^z)^ 



Since, from (j2.1ip 



2 l + |Vu|2 4 l + |Vn[2 2 (l + |Vn[2)2 

+ /?(p- l)Vz.Vn + ez. 



(1 + |Vn|2)2 



we derive from (j2.10p 

(p-2) |Vz|2 ^ ^ (Vz.Vu)2 _ 

which yields, by Young's inequality and the fact that z = |Vn|2, 

- A,z -ip- 2)^^^ + Coz' + 2(m + e)z<C,^-^ + (2.12) 
1 + I Vn|^ I + z 

for some positive constants Cj {j = 0, 1, 2), eventually depending on (3, with the constant 
C2 also depending on \\€u ||oo' Next we introduce the operator J\, defined by 



D^zVu.Vu 

1 + |Vu|2 ' 



A{z) = -/^gZ-{p-2) ^ , , (2.13) 
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which can be written, in local coordinates, as 



A{z) 



kl 



1 + 2^(7 



kl 



E 



ijk 



-fk Im 



nk9'' + g'.{ ) + ip- 2)r^, ^ ^ 



V 



Im 



(2.14) 



where the ajj are uniformly elliptic and bounded and the bi are bounded: indeed, it holds 

min(p - 1, l)g'-^iiij < aijCi^j < max(l,p - l)^*^^^^ • 

Therefore from (j2.12p z is a positive subsolution of an equation of the type 

A{z) + h{z)+g{z)\Vz\^ = f, (2.15) 

where g{z) = — Ci(l + z)"^, h{z) = 2(m + e)z + Cqz^ and / = C2. Since m > 0, g and 
h are increasing functions of the nonnegative variable z, it follows that the comparison 
principle holds between super and sub-solutions of 



AgZ-{p-2)- 



D'^zVu.Vu 



+ Cqz"^ + 2{m + e)z - Ci 



l + z 



Co. 



(2.16) 



Standard computations show that, if A and fj, are positive constants large enough, the 
function 

is a supersolution of (|2.16p . which in addition blows up on dS. We conclude that any 
bounded subsolution of (|2.16|) satisfies z{x) < z{x), and therefore any subsolution by 
replacing S hy {x & S : p{x) > h} and p{x) by p{x) — h. 

Finally, we proved that any u S C'^{S) which is solution of ()2.8p satisfies 



\Vu{x)\ < 



p{x] 



+ Li Vx G S, 



(2.17) 



for some constants Lq, Li depending on ||eti~||oo- Moreover, Lq and Li can be chosen 
uniformly bounded with respect to /3, provided f3 remains in a compact subset of (0, 00). 

To conclude with the estimates on solutions of (12. Sp . it is classical from the theory of 
quasilinear elliptic equations (see e.g. [14j) that local Lipschitz estimates imply local C^'" 
estimates since the equation is smooth and uniformly elliptic. 
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Step 3: Existence for the approximate equation. As in [T7], we consider, for n G N the 
solution Vn,e '■= V of 

D^vVv Vv 

-AgV-ip-2) ^^|^^|2 + Pip - l)\Vv\^ + ev = inS ^^.is) 

v{x) = n on dS, 

By previous steps, the following estimates hold in S. 

0<vn,eix) < -^lnpix)-Mopix) + ^, (2.19) 

|V^;„,(x)| < + (2.20) 

Moreover the sequence {vn,e} is bounded in Cf^'^{S), which ensures the local compactness 
of the gradients. Since n i— > Vn,e is increasing, there exists = lim^^oo Vn,e and is a 
solution of (p:2|) which satisfies ([2J9]) and ([2:20]) . 

Step 4- The ergodic limit. From Step 1, by comparison with Uh and Uf^ defined in 
(j2.6p - (j2.7p (and letting /i — s- 0), we know that their holds in S: 

- ilnp(x) + Mop(x) - ^ <v,{x) < -^In p{x) - Mop{x) + (2.21) 

Therefore ev^ is locally bounded in S. Since Vv^ is locally bounded too in S, €nVe„ 
converges to some constant Aq > for some sequence {e^} in the Qoc-topology of S. We 
fix xo G 5* and set We := Ve{x) — Ve{xo). Because We is locally bounded in Clg^{S) and We 
satisfies 

- Agw, -{p- 2) ^'7;^^"^^^^ + f3{p - l)|Vu;,p + ew, = -ev,{x^) in S (2.22) 

1 + I VWel^ 

the regularity theory for elliptic equations implies that is locally bounded in C^'"(5'). 
Up to an extraction of subsequence, there exists wq = lim„^oo ■, and W{j is a solution of 

-A,z.o-(p-2)^^!!^J^^+/3(p-l)lV^oP = -Ao in 5. (2.23) 

The only question which remains to be proved is that blows-up at the boundary. We 
set ^ 

±{x) = --lnp(x) +Mop(x), 

and get, with same computations as in (j2.4p . 

-Agl -ip- 2)—f-^-f + Pip - 1) I V^|2 + e± 

- . (2-24) 
= I - ^pHp) - 2(P - l)Mo|Vpp) + #(x), 
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where ipfs is a bounded function (depending on /3, Mo). Noticing that |V/)| = 1 in a 
neighborhood of dS, and that eVe{xo) is uniformly bounded, we can choose Mq, po such 
that the function ^ is a subsolution of (j2.22p in {j; G 5 : < p{x) < Po}- Since, whenever 
p{x) = pq, we have Wi:{x) > — cq for some cq > (due to the gradient estimate for v^), and 
since tp — c is stiU a subsolution for any positive constant c, we derive 



yJeix) > — — ln/9(x) + Mop(x) — c Vx s.t. /9(x) < Po- 



(2.25) 



Letting e tend to implies that linix^ds wo{x) = oo. 

Step 5: Uniqueness of the ergodic limit. We claim that there exists a unique constant 
Aq > such that there exists vo £ C'^{S) solution of 



-AgVo-{p-2) — 



D^voVvo-Vvo 
+ |Vt;o|2 



+ P{p-l)\Vvo\ 



-Ao in S 



lim 



x^dS 



vo{x) = oo. 



(2.26) 



To this purpose, it will be useful the following 



Lemma 2.2 A function vo € C'^{S) is solution of \2.26\) if and only if the function 
ujo = e-^^o G C^S) n C{S) is a solution of 



j -divg {{(5^ujI + \yuoo\'^Yl^-^VuJo) = I3\o{P^ojI + |Vwo|^)^/^"^^o in S 
y Wo = on dS. 

Moreover, ujo G C^''^{S) for some 7 > 0, and d^ujo < on dS. 



(2.27) 



Proof. Let vo G C^{S) be a solution of i^Mij- As in the previous steps, considering 
the functions 



^(x) = -4 In p(x) + Mop{x) - M* and 



^\np{x)-Mop{x) + M\ 



which appear to be respectively a sub and a super-solution for (j2.26|) in {x : p[x) < 6} 
for some 6 > small enough (where M* depends on the value of vo on the set {x G S" : 
p{x) = 5}), we obtain 

, , Inp(x) 
vo{x) + 



By the gradient estimates of Step 2, there holds 



< M*. 



|V^;o(x)| < - 



p{x] 



+ Li. 



(2.28) 



(2.29) 



Now set UJO = e-^^o, then coo G T^^'°°(5) nC(5) solves the problem ([227]). By the regular- 
ity theory for degenerate equations of p~Laplacian type (see the Appendix, Theorem lA.ll 
and related references), we can deduce that ujq G C^''^{S). Moreover, since (j2.28p implies 



-f3M* 



< 



^0 

p{x) 



< e 



IBM* 



(2.30) 
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we deduce that d^oJo < — e~^*^ < on dS. As a consequence, since ujq G C^{S) and 
is positive in S, we deduce that problem (j2.27p is uniformly elliptic, so that the classical 
regularity theory applies to give loq G C^'"(5). 

Of course, the converse is also true: given a solution ojo of (|2.27p . clearly vq = — In 
is a solution of ()2.26p . h 

Assume now that there exist two ergodic constants, Ai and A2, associated with two 
solutions vi, V2, and let correspondingly = e~^"' be solutions of (j2.27p . Notice that 
multiplying ()2.27p by loq and integrating on S, we get actually Aq > 0. Thus Aj > and, 
say, A2 > Ai. 

Since u)i/u2 G L°°{S) (from estimate (I2.30p ). we denote 

a = sup — . 

s ^2 

Because equation (|2.27p is homogeneous we can assume that 9 = 1 and either there exists 
xq G S such that uJi{xo) = uj2{xo), Va;i(xo) = Vuj2{xo) and uJi{x) < u)2{x) for x £ S, or 
a;i(x) < a;2(x) for x G 5 and there exists xq G dS such that OuU^i^xq) = dui02{xo). In 
the first case, it turns out that the function z = vi — V2 is nonnegative in S, achieves a 
minimum at xq G S" and satisfies 

A / ^ / ^^ D^z{xo)Vviixo).Vvi{xo) 
-A,z(xo) -{p- 2) ,^|^^^(^^)|2 = A2 - Ai > 0, 

which is impossible because of ellipticity. In the second case, we have du{uJi —^2)ixo) = 0, 
whereas wi — uj2 is negative in 5 and {loi — uj2){xo) = 0. Since the problem ()2.27p is 
uniformly elliptic (recall that the functions tUi satisfy {0^wf + |Vu;jp) > on 5") this 
contradicts Hopf maximum principle. Therefore uji = 0J2, which implies Ai = A2 by the 
equation. Thus the ergodic constant is unique. 

In a similar way one can prove that loq is unique up to a multiplicative constant, and 
so Vq is unique up to an additive constant (as a consequence, the whole sequence w^, 
constructed in Step 4, converges to as e — > 0). 

However, the uniqueness of vq can be proved with a more general argument, concerning 
directly problem (12.26P , which is a variant as well as a generalization of previous uniqueness 
results for explosive solutions. Since it can have its own interest, we present it here. 

First of all, we recall that any function vq solution of (j2.26p satisfies ()2.28p and 
(12:29]) . Moreover, by Lemma Owe have that lvq = e'^^^ G C^(5) and d^ujQ < on dS, 
hence, using that Vvq = —^^S/luq and the estimate ()2.28p we conclude that there exists 
a constant o" > such that, in a neighborhood of dS 

\Vvo\>^. (2.31) 

In addition, it is possible to deduce from (j2.28p - (|2.29p that there exists a constant Cq > 
such that 

\D^vo\<^^ Vxg5. (2.32) 
p'^lx) 
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Indeed, take xq € and let pQ = , where we recall that p{xq) = dist(xo,55). Then 
consider (in a local neighborhood of xq) the rescaled function 

In PQ 



uoiC) = vo{xo + PoO + 



/3 



for C G -6(0,1). Note that p{xo + PoO ^ (po,3po) so that (l2:29|) and (123X1) imply | < 
I-DuqI ^ Lq+Li pq. Since f;o is a solution of (j2.26|) . a simple scaling in the local coordinates 
gives that uq is a solution of 

-A,no - (p - 2):^^^^^^^ + P{p - l)|VnoP = -Ao pI for ? G B{0, 1) 

with a slight abuse of notation since now, in local coordinates, the derivatives are taken 
with respect to the variable ^. Since the second order operator is uniformly elliptic (as 
in (j2.13p - (j2.14p ). by the classical regularity theory (e.g. see [H], Theorem 13.6 to deduce 
the Holder estimates for Duq and then apply the Schauder estimates. Chapter 6) we have 
that 

where C is a constant depending on sup^^o.i) (l^ol + l^'i^ol)- Using the estimates (j2.28P " 
(j2.29p we can bound this last quantity only depending on M*, Lq, Li, hence we conclude 
that \D'^uo{0)\ < C, which gives (|232D . 

Now, take two solutions vi, V2 of (j2.26p corresponding to Ai, A2 with, say, Ai < A2. 
We adapt now an argument in [17J: consider the function v = 0v2, for 9 < 1, and compute 



D^vVv.Vv ^, ..,„.|2 

~ — uA 



+(1 - 9')9{p - 2) (T^^^ffg^ll^ - (1 - 9)epip - 1)|V.2 



2 

2 



Using (|2:29]) . (|232]) and dOT]) . we know that 

D'^V2Vv2.'Vv2 



{1 + \VV2\^){1+9^\VV2\^) 



<C in 5 



hence v satisfies 



-A,f;-(p-2) ^'^y^.y^% /?(p-l)|Vl)[^< 

Thanks to (12. Sip , we deduce that there exists 5 > 0, independent on 9, such that v satisfies 

-A,v -ip- 2)^^^^ + (3ip - l)|Vi)p < -Ai 
1 + I vu|^ 

in {x G S : p{x) < 6}. However, from the estimate (|2.28p which holds for vi and V2 we 
have that vi — v ^ +00 as p{x) — 0, hence vi — v has a minimum in {x G S : p{x) < 6} 
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and, by standard maximum principle, it is reached when p{x) = 5. Letting 9 ^ 1, we 
conclude that 

min{(t>i — V2){x), x : p{x) < (5} = min{(t;i — V2){x), x : p{x) = 5} . 

On the other hand, looking at the equations of f i, ^2 in {x G 5 : p{x) > (5}, we also know 
(again by maximum principle) that 

min{(fi — V2){x), x : p{x) > 5} = min{(i)i — U2)(x), x : p{x) = 5} 

hence vi — V2 should have a global minimum reached at a point xq € 5* such that p{xq) = 5. 
Since xq lies inside the domain, and the function z = vi — V2 satisfies a smooth elliptic 
equation around xq, using the strong maximum principle we conclude that vi — V2 is 
constant. This proves the uniqueness, up to a constant, of the solution of (j2.26p . and at 
the same time also the uniqueness of the ergodic constant (Ai = A2, as we already proved 
before) . □ 

Remark 2.3 The argument used in the last step of the previous proof also provides a 
general uniqueness result for explosive solutions of 

-AgV-{p-2)^^^^^ + P{p-l)\Vv\^ + ev = f inS ^2.33) 

lim^^QSv{x) = 00. 

Precisely, if / is a Lipschitz function, and e > 0, the problem (|2.33|) has a unique solution 
V G C^(S'). To our knowledge, such a result is new even in the euclidean setting M = R^. 



We proceed now studying how the ergodic constant depends on /?, which will lead 
to the proof of Theorem B. 

Proposition 2.4 Under the assumptions of Theorem \2.1{ the mapping (3 ^ \p is contin- 
uous and decreasing from (0, 00) in (0,oo), and it verifies 

lim Afl = 00. (2.34) 

Proof. Step 1: the monotonicity. Let < /3i < /32 and let v^^i and Ve,2 be the corresponding 
solutions of (12. 2[) with /? respectively replaced by j3i and /?2. Since the Ve,i are limit 
of solutions with finite boundary value there holds v^^i > 'Ue,2 by comparison principle. 
Therefore 

A/9, := lim ev^ 1 > Aa, := lim ev^ 2- 

Next, if we assume that there exist Pi {i = 1, 2) such that < /3i < /32 and = = A 
and if uji and uj2 are the corresponding solutions of ()2.27p with f3 = (3i and A = A/j^ = Xjs^, 
then ([2:28|) implies 

m~^p{x) < uJi < mp{x) Vx € S, 
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for some m > 0. Set Co = uj^'^^^^ , then 

-divg [{PICj^ + |VlD|2)p/2-1VlD) - hKPloj'^ + |VcD|2)P/2-i 



UJ 



P2\ f P2Y ^ .(p-l){/32//3i-l) (^^2^2 ^ |Vwi|2)(P-2)/2 |V^lP _ (2-35) 



Therefore cj is a strict sub-solution. By homogeneity, and since di^oj vanishes on dS, we 
can assume that oj < uj2, that there exists xq G S such that u}{xo) = u}2{xo) and the 
coincidence set of uj and UI2 is a subset of 5. Let 

z = — — (lna;2 — Inu;) = f2 — v. 

Then z < 0, it is not identically zero, z{xq) = and z{x) — > —00 as p{x) dS. Because 

D'^V2Vv2.Vv2 



-AgV2-{p-2y- 



1 + iVval^ 



+ p2ip - 1)\^V2\' 



< -AgV -{p- 2)———^ + P2{P - l)|Vii|2 



'^^ l + |Vi)|2 

developing this inequality, we obtain that, at x = xq, there holds 

D'^Z\IV2SIV2 



+(P - 2) 



-^^^-(^-')-+|V.2P 
D^vVv.Vv D^vVv2.Vv2' 



+ f32{p-l) [|V?;2p - |VS|2] < 



l + |Vf}|2 i + lVi-aP 

Since ^2 are in S", the strong maximum principle yields a contradiction. Therefore 
/? I— > is decreasing. 

S'tep 2: the continuity. Let {/?«} be a positive sequence such that /3„ — > /3o and Vf^^ be the 
corresponding solution of 

'""-SW^^""'^-"'^"-'^-^- '"^ (2.30, 

lim^._,a5 (x) = 00, 

and let v^^js^ be the corresponding solutions of (|2.2p with /3 = Since efe^/3„ remains 
locally bounded in S when /?„ remains in a compact subset of (0, cxd) and converges to Xp^ 
locally uniformly as e ^ 0, the set {Xp^} is bounded. Up to a subsequence (not relabeled) 
we can assume that A^^ ^ A as n — > 00. Thanks to (j2.28p and (j2.29p . there holds 



^/3„ + 



In p{x) 



f3n 



< Co and {"^vpj < 



Ci 

p{x)' 



(2.37) 



for some constants Cq, Ci, hence the sequence {v/3„} remains locally bounded in Wl^^{S) 



and, therfore, in Cf^'^{S). Up to a subsequence Vf^^ — > z; in C2^^(S'), and u is a solution of 

D'^vVv.Vv 



+ /3o(p-l)|V7)| 



-A in S 



lima;^a5i;(x) = 00. 
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By uniqueness of the ergodic limit, A = Xp^, and Ayg,^ Xp^ for the whole sequence. 
Step 3: {2. 34^ holds. Let be a positive solution of 

r -divg {{P^Lo^ + |Vcj|2)p/2-1Vl^) = ^A/3(/?V + |Vc^[2)P/2-i,^ in S 
< (2.38) 

w = on 95. 



We normalize to by 



\VLo\Pdvg = 1. 



Therefore, if /U^ is the first eigenvalue of —diVgdV . |^~^V .) in Wq'^{S), there holds 

/ \ojFdv„ < — . 
J /^s 

Multiplying (I2.38P by uj and integrating over 5 yields to 

[ ip^uj^ + \Vuj\y/^dvg = (5{Xp + I3)f {(3^uj^ + \Vuj\^f/'^-^uj^dvg. (2.39) 
Js Js 

Clearly 

/ iP^UJ^ + |Vwl2)P/2d-t;g > [ \Vuj\PdVg = 1. 

Js Js 

If p > 2, 

/ (/32c^2 ^ |Vw|2)P/2-1^2^^^ < 2P/2-2 /" (ujP + UJ^\Vuj\P-^)dVg 

Js Js 

< 2P/2-2 (l + ^) j UjPdVg + 2P/2-2 (l - ^) j \Vuj\P)dVg 



This implies 



If 1< p < 2, 



Therefore 



< C „ 

— p,S 



/3(A/3 + /3) > ^ ^ > -1- - p. (2.40) 



^5 



/J^-^A^ + P)>fis^^P>^-P- (2-41) 
Clearly (|2:i0]) and (l2:iT]l imply ([23^ . □ 

Remark. Using the uniform ellipticity and the maximum principle, (j2.40p and (j2.4ip can 

C 

possibly be improved in A^ > — . 
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We have now all the ingredients for the proof of Theorem B. 

Proof of Theorem B. If we set lo = e~^^ where v is the solution of (j2.ip . then 
Lv is defined up to a multiplicative constant and satisfies (j2.38p . By Lemma 12.21 oj G 
C^{S) n C^(S'). Therefore the Theorem is obtained if we can prove that there exists a 
unique (3 := j3g > such that 

\p = f3{p-l)+p-d-l. (2.42) 
But the mapping /? i— > — f3{p — 1) is continuous and decreasing on (0, oo). Clearly 

lim A/3 — /3(p — 1) = — oo, 

and 

lim Afl — /3(p — 1) = cxD, 
by Proposition 12. 4[ The results follows by continuity. □ 

3 The regular case and Tolksdorf 's result 

If /? < 0, the equation satisfied by a a separable p-harmonic function u under the form 
p.ip is unchanged. However, if we set j3 = — then ()1.4p turns into 

- div [{P'^uj + I V'cj|2)P/2-iv'a;) = P{p(p -l)+N- p)CP'^uj + \Vu\y'^~^u. (3.1) 

Furthermore, if a solution oj of (13. ip in 5" C S^^^ exists which vanishes on dS, then 
f3{p — 1)+A^— p>Oby multiplying by uo and integration over S. By setting 

In a; 



then V satisfies 

' -div ((1 + \Vv\Y^~^ V'?;) + [3{p - 1) (1 + \Vv\^Y'^~^ \Vv\^ 
< = -0{p - 1) + N - p) {1 + \Vv\'^)'^^^''^ inS (3.2) 

lim^^QS v{a) = oo. 

In the general setting of a Riemannian manifold. Theorem 12.11 and Proposition 12.41 are 
valid with /? replaced by p. The proof of Theorem B holds except that (j2.42p is replaced 
by 

X^ = (3ip-l) + d+l-p. (3.3) 

Because the function /3 i-^ A^ — P{p — 1) is unchanged, the proof of Theorem B applies and 

shows that there exists a unique P := (3^ > such that (j3.3p holds. Consequently we have 
proved the following result which contains Tolksdorf 's initial result if {M,g) = {S^~^ , go). 



17 



Corollary 3.1 Under the assumptions of Theorem \2.1\ there exists a unique (3 := (3g > 
such that the problem 

J -divg ^ I Vw[2)p/2-1Vl^) = /3 - 1) + d + 1 - (^2^2 ^ I Vu;|2)p/2-i^ in S 

[ w = on 95, 

(3.4) 

admits a positive solution uj E C'^{S) fl C^(S'). Furthermore oj is unique up to an homoth- 
ethy. 

A Appendix 

We prove here the C^''^ regularity up to the boundary, stated in Lemma [221 for solutions 
of degenerate equations in divergence form 

{—div(a(x,u,'Vu))=B(x,u,Vu) in S 
(A.l) 
u = on dS. 

We will assume that a{x,s,(,) satisfies the following conditions: there exist constants A, 
A, /? > 0, and a G (0, 1], and a continuous function fi : S x M ^ M such that, for every 
s, i G M, for every ^, r/ € M^, and a.e. x G 

?^{x,s,Omrij>X{fi{x,sf + \C\^)'^\v\\ (A.2) 



<A(Ma^,s)2 + |eP)^, (A.3) 

ia(x, - a(y, t, 01 < /? (1 + ler ' + ler [ix - yr + Is - ri, (a.4) 

\B{x,s,0\<P{l + m- (A.5) 
The model we have in mind is clearly 



P-2, 



a{x,u,Vu) = {n{x,u) + |Vu| ) 2 

where p > 1, and the function n{x,s) is Lipschitz (or possibly Holder) continuous. In 
many cases, as in the proof of Lemma 12.2^ the a priori information that u is Lipschitz (or 
Holder) continuous could allow us to consider only the case jj, = n{x). 

The C'^'^ estimates, or similar kind of regularity results, are by now classical since 
the works of E. DiBenedetto [9] and P. Tolksdorf [20] for the p-Laplace equation: as far 
as the global regularity, up to the boundary, is concerned, we refer to the works of G. 
Lieberman (e.g. [18]) or to [8]. Despite a large amount of literature available, it seems 
that no exact reference applies to our model, so that, for the sake of completeness, we feel 
like giving a proof of this result, at least detailing the possible slight modifications in order 
that previous results can be generalized. To this purpose, we observe that while the case 
p > 2 is somehow contained, if not in previous statements, at least in previous arguments 
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(specifically, we refer to [18]), this seems not sure for the case p < 2 because of our growth 
assumption ()A.4p (roughly speaking, the (x, s)-derivatives may grow like |'^|''~^)- Finally, 
we note that the next result would still hold for a nonhomogeneous boundary condition 
(u = (/3 on dS) provided ip belongs to C^'" {dS). 

Theorem A.l Let S be a bounded C^'" domain in M^, and assume that M.^j -i Pl.^)) 
hold true. If u is a bounded weak solution of (HHP, then there exists 7 G (0, 1) such that 
u € C^''^{S) and moreover 

W^Wc^riis) <C{k/X,a,\\u\\oo,p,N,S) . 

Proof. Because our specific interest is in the boundary estimate, we only prove the 
regularity of u around a point xq € dS (the inner regularity is treated in the same manner). 
Up to straightening the boundary, we can assume that locally dS = {x : xjy = 0} and 
S = {x : XN > 0}. 

We follow the standard approach via perturbation argument. We denote Bji = {x : 
\x — xo\ < R}, = Bji n S, and consider the solution v of 

{—div {a{xQ, u{xo), Vv)) =0 in B^ 
V = u on dB^. 

Problem ()A.6P has a unique solution v G W^'^{B~^). Due to assumptions ()A.2p - ()A.3p . the 
estimates concerning v are well-established ([9], [20], [IE]). In particular, from Lemma 5 
in [18] we have, for some cr > 0, 

r^-/_^ f ,„ \' R 



oscVv <C [-] R-'^ / \Vv\Pdx Vr < - (A.7) 

B+ \RJ y Jb+ J 2 

where C, here and after, depends only on the constants appearing in the hypotheses and 
possibly on |j^||oo5 i^^ particular through the quantity sup{|/i(x, s)| ,x G /S, \s\ ^ ||t^||oo}' 
Moreover, since a(x,s,^)^ > c{\^\p — \fi\^), one easily deduces from (1A.6[1 . using v — u as 
test function and Young's inequality, that 

\Vv\P dx < C (^1 + J ^\Vu\Pdx^ . (A.8) 

Finally, the maximum principle gives inf u < v < sup n, which yields 



R 



osc V < osc u . (A. 9) 

Br 

Now take u — v as test function both in (lA.ip (restricted to B^) and in (IA.6P to obtain 

/ a{x,u,Vu)V {u — v) dx— I a{xo,u{xo),'Vv)'V{u — v)dx= / B{x,u,Vu){u — v) dx . 
J b\ J si J Bi 
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Denote Dy : = {x e : |Vii| < \Vv\} and Du : = {x £ : |V-u| < |Vn|}: hence we 
have 

[a{x, u, Vu) — a{x, u, Vv)] V(n — v) dx 
+ [a{xo,u{xo),Vu) — a(xo, u{xo), Vv)] V{u — v) dx 

(A. 10) 

= [a{xQ,u{xQ)^Vv) — a{x,u,Vv)]V {u — v) dx 
+ [a{xQ,u{xo),Vu) — a{x,u,Vu)]S/{u — v) dx + f^+B{x,u,Vu){u — v)dx 

Using ()A.4j) and the definition of Dy, we have 

Jj^ [a{xo,u{xQ), Vf ) — a(x, u, Vv)] V{u — v) dx 

< 213 (1 + \Vv\P~'^ + \Vv\P-'^) \Vv\ [\x - xoT + \u{x) - u(xo)|"] dx 

< C[R^ + (oscu)"] /^^ (1 + \Vv\P) dx 

B+ 

Similarly we estimate the second term in the right hand side of (jA.lOp . and using also 
()A.5P we deduce 

[a{x, u, Vu) — a{x, u, Vv)] V(n — v) dx 
+ Jjj [a(xo, u{xo),Vu) — a(xo, u{xo), Vv)] V{u — v) dx 

< C[R'' + (osc n)" + osc u] L+ (1 + \Vv\p + \Vu\p) dx , 

where we used that osc(ti — v) <2 osc u thanks to (|A.9p . 

Br b+ 

Now, in both terms in the left hand side we use (|A.2p which implies, for every (x, s, ^): 

[a(x, s, - a{x, s, t?)] c(A) (^(x, s)' + [Cp + |r?P) \^ - ■ (A.ll) 

If p < 2 we get (recall that the generic constant C may depend on ||n||oo) 
[a{x,u,Vu) — a{x,u,Vv)]V{u — v) dx 
+ Jj^ [^(xo, u{xo), Vu) — a(xo, u{xo), Vv)] V(n — v) dx 

> CJd^ud^ [1 + |Vn|2 + |V^;|2] |V(n - v)\' dx 
hence using Holder inequality we end up with 

\V{u-v)\Pdx < C[i?" + (oscn)" + oscu]9 {1 + \Vv\p + \Vu\p) dx , 

^ Br B^ ^ 

with Q' = |. If p > 2 we simply get rid of the term ^u^ in (jA.lip and obtain the same 
inequality with g = 1. Therefore, using also (TO]) , we conclude for any p > \ 

j \V {u -v)\Pdx<C[R°' + {osc u)°' + osc uf j^+{l + \Vu\P) dx ^^^2) 

^ Br B^ ^ 
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with q = min(l, |). 

Starting from the inequahty ()A.12p it is possible to deduce the Holder regularity of 
Vti following well-known arguments. In particular, if u is Lipschitz continuous (as in our 
application in Lemma l2.2p the conclusion is straightforward, since ()A.12p implies 



Jb+ 



and (lA77ll - (IX8l) give that oscVu < C {r / Rf . 

Then, defining (F),^ = F{y)dy for F = Vu or Vv^ we deduce 

|Vn - {Vu)X dx<C [/^+ \Vn - Vvf + \Vv - {Vv^ 

and choosing R = for some suitable 9 < 1 the conclusion follows from the results of 
Campanato [7J. 

In the general case, i.e. when a Lipschitz estimate on u is not available, one need 
further work to estimate the right hand side of ()A.12p . For this purpose, starting from 
(IA.12h . we can follow the arguments of G. Lieberman ([18j, Section 3) and still get at the 
conclusion. ■ 
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